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We propose an implementation of quantum logic gates via 
virtual vibrational excitations in an ion trap quantum com- 
puter. Transition paths involving unpopulated, vibrational 
states interfere destructively to eliminate the dependence of 
rates and revolution frequencies on vibrational quantum num- 
bers. As a consequence quantum computation becomes feasi- 
ble with ions whos vibrations are strongly coupled to a ther- 
mal reservoir. 
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Recently, methods to entangle states of two or several 
quantum systems in controlled ways have become subject 
of intense studies. Such methods may find applications 
in fundamental tests of quantum physics JjJ and in preci- 
sion spectroscopy 0]; and they offer fundamentally new 
possibilities in quantum communication and computing 
||. A major obstacle to these efforts is decoherence of 
the relevant quantum states. In many proposed imple- 
mentations of quantum computation the quantum bits 
(qubits) are stored in physical degrees of freedom with 
long coherence times, like nuclear spins, and decoherence 
is primarily due to the environment interacting with the 
channel used to perform logic gates between qubits jjj. 
In this Letter we present a scheme which is insensitive to 
the interaction between the quantum channel and the en- 
vironment. Specifically, we consider an implementation 
of quantum computation in an ion trap, but we hope 
to stimulate similar ideas to reduce decoherence in other 
physical implementations. 

The ion trap was originally proposed by Cirac and 
Zoller H as a system with good experimental (optical) 
access and control of the quantum degrees of freedom, 
and they suggested an implementation of the necessary 
ingredients in terms of one-bit and two-bit operations to 
carry out quantum computation. In the ion trap com- 
puter, qubits are represented by internal states of the 
ions. The number of qubits equals the number of ions, 
and this system is scalable to the problem size in contrast 
to NMR quantum computation which is only applicable 
with a limited number of qubits . 

The ion trap method || uses collective spatial vibra- 
tions for communication between ions, and it requires 
that the system is restricted to the joint motional ground 
state of the ions. For two ions this has recently been 
accomplished H]. We present an alternative implemen- 
tation of quantum gates that is both insensitive to the 
vibrational state and robust against changes in the vi- 



brational motion occurring during operation, as long as 
the ions are in the Lamb-Dicke regime, i.e., their spatial 
excursions are restricted to a small fraction of the wave- 
length of the exciting radiation. Our mechanism relies 
on features of quantum mechanics that are often respon- 
sible for "paradoxical effects" : i) The vibrational degrees 
of freedom used for communication in our scheme only 
enter virtually i.e., although they are crucial as interme- 
diate states in our processes, we never transfer population 
to states with different vibrational excitation, ii) Tran- 
sition paths involving different unpopulated, vibrational 
states interfere destructively to eliminate the dependence 
of rates and revolution frequencies on vibrational quan- 
tum numbers. 

Like in the original ion trap scheme ||, we address 
each ion with a single laser, but quantum logic gates in- 
volving two ions are performed through off-resonant laser 
pulses. For the laser addressing the first ion, we choose a 
detuning close to the upper sideband, i.e., close to reso- 
nance with a joint vibrational and internal excitation of 
the ion. We choose the detuning of the laser addressing 
the second ion to be the negative of the detuning of the 
first laser, see Fig. [I] (a). This laser setting couples the 
states \ggn) <-> {|egn+ 1), \gen — 1)} <-> |een), where the 
first (second) letter denotes the internal state e or g of 
the first (second) ion and n is the quantum number for 
the relevant vibrational mode of the trap. We choose the 
detuning from the sideband so large that the intermedi- 
ate states \egn+ 1) and \gen— 1) are not populated in the 
process. As we shall show below, the internal state tran- 
sition is insensitive to the vibrational quantum number 
n, and it may be applied even with ions which exchange 
vibrational energy with a surrounding reservoir. 

If we tune the lasers sufficiently close to the sidebands 
we can neglect all other vibrational modes and concen- 
trate on one collective degree of vibrational excitation of 
the ions ||. In this case our system can be described by 
the following Hamiltonian 



H — H n + H int 
H = Tiv(a)a + 1/2) + %uj. 
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where v is the frequency and a' and a are the ladder op- 
erators of the vibrational mode and Tiw eg is the energy 
difference between the internal states e and g. Pauli ma- 
trices (Ti represent the internal degrees of freedom for the 
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i'th ion, and u}% and Vli are the frequency and Rabi fre- 
quency of the laser addressing the i'th ion. [ In a practical 
realization, one might use Raman transitions between low 
lying states of the ions due to their long coherence time. 
By appropriate redefinition of the symbols our formalism 
also describes this implementation 0.] We consider an 
ion trap operating in the Lamb-Dicke limit, i.e. the ions 
are cooled to a regime with vibrational quantum num- 
bers n ensuring that r\i\Jn + 1 is well below unity (Note 
that this may still allow n- values well above unity.) In 
our analytical calculations we use an expansion of H mt 
to second order in rji , but in our numerical treatment we 
apply the exact Hamiltonian (|l|). 

We wish to perform an operation on the mutual state of 
two ions 1 and 2 selected freely within the string of ions, 
and we assume that 771 = f]i = r\ and f^i = = n. With 
the choice of detunings described above, the only energy 
conserving transitions are between \ggn) and \een). The 
Rabi frequency fl for the transition between these states, 
via intermediate states m, can be determined in second 
order perturbation theory, 



( een \ H mt | to) (to | H int | ggn) 
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where the laser energy TiUi is the energy of the laser ad- 
dressing the ion which is excited in the intermediate state 
|m). If we restrict the sum to \egn+ 1) and \gen— 1), we 
get 



= - 



m 2 



(3) 



where § = u>i — uj eg is the detuning of the laser addressing 
the first ion pQj . 

The remarkable feature in Eq. ^ is that it contains 
no dependence on the vibrational quantum number n. 
This is due to interference between the two paths indi- 
cated in Fig. |l|(b). If we take the path where ion No. 
1 is excited first, we have a factor of n + 1 appearing in 
the enumerator (^/n + 1 from raising and \Jn + 1 from 
lowering the vibrational quantum number). In the other 
path we obtain a factor of n. Due to the opposite de- 
tunings, the denominators in Eq. (||) have opposite signs 
and the n dependence disappears when the two terms 
are subtracted. The coherent evolution of the internal 
atomic state is thus insensitive to the vibrational quan- 
tum numbers, and it may be observed with ions in any 
superposition or mixture of vibrational states. 

From the above arguments we expect to see perfect 
sinusoidal oscillations between the population of the in- 
ternal states I gg) and |ee). To confirm the validity of 
our pertubative analysis we have performed a direct nu- 
merical integration of the Scrodinger equation with the 
Hamiltonian ([!]) to all orders in to We have considered 
a situation, where both ions are initially in the internal 



ground state. For the vibrational state, we have investi- 
gated a number of different states, including Fock, coher- 
ent and thermal states, all yielding qualitatively similar 
results. The outcome of the computation for a coherent 
state of vibrational motion can be seen in Fig. |^, where 
we show the evolution of relevant terms of the atomic in- 
ternal state density matrix pij t kt — Tr n (p\kl)(ij\), where 
2, j, k, I = e or g, and where Tr n denotes the partial trace 
over the unobserved vibrational degrees of freedom. The 
figure clearly shows that we have Rabi oscillations be- 
tween the atomic states \gg) and |ee), and the values 
of the off diagonal element p gg ,ee confirm that we have 
a coherent evolution of the internal atomic state which 
is not entangled with the vibrational motion. Superim- 
posed on the sinusoidal curves are small oscillations with 
a high frequency due to off resonant couplings of the type 
\ggn) -> \egn + l),\ggn) \gen - 1) and \ggn) -> \egn). 
The magnitude of these oscillation and the deviation from 
ideal transfer between \gg) and |ee) can be suppressed by 
decreasing n. 

The analysis given so far is sufficient for creation of 
internal state entanglement, completely decoupled from 
the external motion of the ions. By optical pumping we 
can prepare a state p = \gg)(gg\ <8> p V ib, and if we apply 
radiation fields corresponding to a pulse of duration T = 
^j^j- , our system is described by the density operator p — 

<8> Pvib, where is & maximally entangled EPR- 
statc -^Qgg) - i\ee)). 

Since the states \eg) and \ge) do not fulfill any res- 
onance condition one might expect that they are unaf- 
fected by the laser pulses. Due to n-dependent perturba- 
tions of the energy levels by the lasers this is, however, 
not the case. Keeping only the most important terms, 
we get the energy shifts 
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The energy shifts of the \een) and \ggn) are identical and 
independent of n, but since the energy shifts of \egn) 
and \gen), depend on the vibrational quantum number, 
the time evolution introduces phase factors e - lAE egnt/h ^ 
which depend on n, and e.g., at the time t = Ti nv — 
2 ^iq-P where \gg) end |ee) are inverted, factors of (—1)" 
will tend to extinguish the coherence between internal 
states |ee) and \eg). This coherence can be restored by 
a trick resembling photon echoes [0. Notice that the 
n dependent part of AE egn is minus the n dependent 
part of AE gen . If at any time T/2 we change the sign 
of the laser detuning 5, phase components proportional 
to n will begin to rotate in the opposite direction and at 
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time T we will have a revival of the coherence. This is 
confirmed by our numerical solution of the Schrodinger 
equation presented in Fig. [| where we change the laser 
detunings at the time Tj„„/2 and at the time Ti nv we have 
completed the transfer -±=(\gg) + \eg)) — ► -^(— i|ee) + 

\eg)). 

No particularly demanding assumptions have been 
made for the experimental parameters. With a vibra- 
tional frequency i>/2tt — 200 kHz, the transition shown 
in Fig. ||, require Rabi frequencies Q/2tt of modest 20 
kHz, and the evolution from \gg) to |ee) is accomplished 
in 5 ms. To be relevant for real computational tasks it 
is necessary that our evolution is robust against decoher- 
ence effects on this long time scale. An important source 
of decoherence is heating of the vibrational motion, and 
it is a major asset of our proposal that it can be made 
insensitive to the interaction with the environment: The 
arguments leading to Eq. (||) do not require that the ions 
remain in the same vibrational state, and the coherent 
oscillation from \gg) to |ee) may still be observed when 
the vibrational motion exchange energy with a thermal 
reservoir. The photon echo trick, however, is sensitive to 
heating: If the vibrational quantum number n change its 
value at the time T/2 where the detunings are inverted, 
the second half of the gate, will no longer revert the phase 
evolution due to the new value of n and coherence is lost. 
If instead the detunings are inverted N times during a 
gate, the erroneous phase will only be on the order of 
the phase evolution in time T/N, and the effect of the 
heating is reduced. 

Rather than inverting the detunings a large number of 
times, we suggest to continuously apply lasers with both 
detunings ±5 on both ions. With two fields of opposite 
detunings and identical Rabi frequency Q there are two 
contributing paths in addition to the two paths in Fig. 
|]. The contribution from the two additional paths are 
identical to the two original paths and the only modifi- 
cations to the | gg) <-> |ee) Rabi frequency in Eq. ([|) is 
multiplication by a factor of two. With bichromatic fields 
there also exists a resonant transition from \eg) to \ge). 
The Rabi frequency of this transition is the negative of 
the Rabi frequency from \gg) to |ee) and the evolution 
will be described by 
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To validate that the evolution in Eq. (g) is in fact 
stable against heating we introduce a thermal reservoir 
described by relaxation operators c\ = yT(l + n t i lcrm )a 



and C2 = VrTithermfl i where T characterizes the strength 
of the interaction and rithcrm is the mean vibrational 
number in thermal equilibrium. We analyse the dynam- 
ics of the system using Monte Carlo wavefunctions |l2) , 
which evolve with a non Hermitian Hamiltonian inter- 
rupted by jumps at random times. The result of the com- 
putation can be seen in Fig. [|, where we show (left) the 
result of a single Monte Carlo realization with quantum 
jumps indicated by arrows and (right) the average over 
10 realizations. In the figure we have chosen ntherm = 2. 
This rather low value could represent a heating mecha- 
nism counteracted by laser cooling on a particular ion re- 
served for this purpose. In the simulations 34 vibrational 
quanta are exchanged with the reservoir on average, and 
we wish to emphasize that with the proposed scheme the 
gate is almost unaffected even though the duration of 
the gate is much longer than the coherence time of the 
channel used to communicate between the qubits. 

It has been proven that any unitary evolution involving 
a number of qubits can be constructed using only single 
qubit operations and a simple universal quantum gate 
like for instance the two-qubit control- not operation [jl3| . 
With the evolution described by Eq. (||) a control-not 
operation is created by the following sequence of opera- 
tions: Pi, P 2 ~ 1 , H 2 , R, Pi, Hi, Pi, R and P 2 , where R is 
the evolution in Eq. (||) with T — Pi is a ir/2 phase 

change of |e) in ion i, and Hi is a Hadamard transforma- 
tion on ion i. With ready access to one-qubit operations 
in the ion trap we thus have available the ingredients for 
successful quantum computation. 

We note that with only two ions in the trap, the use of 
bichromatic light leads to the evolution (H) even without 
individual optical access, and if many ions are illumi- 
nated by the same field, multi particle entanglement is 
created As a further remarkable feature of our im- 
plementation, we note how easy it is to simultaneously 
operate on different pairs of ions: If we wish to apply 
our gate on the pairs (i,j) and (k,l), we simply illumi- 
nate ions i and j with fields of detunings ±<5y and ions 
k and I with another pair of detunings iSki- The only 
resonant transitions are the ones of the two desired gates 
and, however mind boggling this may be, the virtually 
excited but never populated vibrational mode has been 
used for two processes at the same time. 

Note added: Since submission of this work, two pro- 
posals for computing with vibrationally excited ions have 
appeared. pl| uses widely separated vibrational states, 
and can hardly be generalized beyond 2 or 3 qubits; |l6| ] 
involves dynamics which is conditioned on the vibrational 
state and is not applicable during heating. 
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FIG. 2. Rabi oscillations between |<?<?) and |ee). The figure 
shows the time evolution of the internal atomic density ma- 
trix elements p 99 , 99 (full line), p ee ,ee (long dashed line) and 
Im(p 99iee ) (short dashed line). The magnitude of Re{p gg ^ e ) 
is below 0.03 and is not shown. In the initial state, the ions 
are in the internal ground state and a coherent vibrational 
state with mean excitation n = 2. Parameters are 8 — 0.90 v, 
Q, = 0.10 v and r) = 0.10 
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FIG. 3. Revival of coherence at Ti nv . Initially ion 1 is 
in a superposition -±=(\g) + |e)) and ion 2 is in \g). The sign 
of 8 is changed at Ti nv /2 (left arrows), to ensure the perfect 
transition at Ti nv (right arrows) . The full lines represent pop- 
ulations of |ee) and \eg). The dotted and the dashed curves 
represent the real and imaginary part of the coherence p e e,e 9 - 
Parameters are 8 — 0.90 v, Q = 0.05 v and 77 = 0.1 
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FIG. 1. Energy levels and laser detunings. (a) Two ions 
with quantised vibrational motion are illuminated with lasers 
detuned close to the upper and lower sidebands, (b) The ions 
oscillate in collective vibrational modes, and two interfering 
transition paths are identified. 
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FIG. 4. Rabi oscillations in a heating trap. The left panel 
shows the result of a single Monte Carlo realization with a 
total of 39 jumps occuring at times indicated by the arrows. 
The right panel is an average over 10 realizations. The curves 
represent p 99 , 99 (full line), p ee , ee (long dashed) and Jm(p 99i<3e ) 
(short dashed). Parameters are S = 0.90 v, Q = 0.10 v, 
r) — 0.1, r = 2 * 10~ 4 v and n tcrm = 2. 
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